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Abstract 

6v 

D | In this paper we discuss massive gravity in de Sitter space via gravita- 

tional Higgs mechanism, which provides a nonlinear definition thereof. The 
Higgs scalars are described by a nonlinear sigma model, which includes higher 
OJ derivative terms required to obtain the Fierz-Pauli mass term. Using the afore- 

iyv ■ said non-perturbative definition, we address appearance of an enhanced local 

symmetry and a null norm state in the linearized massive gravity in de Sitter 
space at the special value of the graviton mass to the Hubble parameter ratio. 
By studying full non-perturbative equations of motion, we argue that there is 
f^ ■ no enhanced symmetry in the full nonlinear theory. We then argue that in the 

^> \ full nonlinear theory no null norm state is expected to arise at the aforesaid 

special value. This suggests that no ghost might be present for lower gravi- 
ton mass values and the full nonlinear theory might be unitary for all values 
of the graviton mass and the Hubble parameter with no van Dam-Veltman- 
Zakharov discontinuity. We argue that this is indeed the case by studying full 
nonlinear Hamiltonian for the relevant conformal and helicity-0 longitudinal 
modes. In particular, we argue that no negative norm state is present in the 
full nonlinear theory. 



(N 



ctf 



1 Email : alberto . iglesias@physik.uni-muenchen.de 

2 Email: zura@quantigic.com 

3 DISCLAIMER: This address is used by the corresponding author for no purpose other than 
to indicate his professional affiliation as is customary in scientific publications. In particular, the 
contents of this paper are limited to Theoretical Physics, have no commercial or other such value, 
are not intended as an investment, legal, tax or any other such advice, and in no way represent 
views of Quantigic® Solutions LLC, the website www.quantigic.com or any of their other affiliates. 



1 Introduction and Summary 

Gauge interactions are mediated by massless particles, gauge vector bosons. Upon 
spontaneous breaking of gauge symmetry, gauge bosons acquire mass via Higgs 
mechanism. Similarly, components of the massless particle associated with gen- 
eral coordinate reparametrization invariance, the graviton, could acquire mass via 
gravitational Higgs mechanism ppj upon spontaneous breaking of diffeomorphism in- 
variance by scalar vacuum expectation valuesj. The gravitational Higgs mechanism 
was revisited in the context of obtaining massive gravity directly in four dimensions 
in|23ll21ll251l2Sll2Zll251l2Sll3QlEIlE2]- A general Lorentz invariant mass term for 
the graviton h MN is of the form 

- ^ [h MN h MN - K{h$f\ , (1) 

where k is a dimensionless parameter. Unitarity requires that the graviton mass 
term be of the Fierz-Pauli form with k = 1 [33J. Massive gravity in Minkowski 
space with k — 1 can be obtained via gravitational Higgs mechanism by including 
higher derivative terms in the scalar sector and appropriately tuning the cosmological 
constant against the higher derivative couplings |26j. 

The framework of [26J provides a ghost-free, nonlinear and fully covariant def- 
inition of massive gravity in Minkowski space via gravitational Higgs mechanism 
with spontaneously (as opposed to explicitly) broken diffeomorphisms. In this pa- 
per we discuss massive gravity in de Sitter space via gravitational Higgs mechanism, 
which provides a nonlinear definition thereof. This is achieved by coupling gravity to 
scalars, whose vacuum expectation values result in spontaneous breaking of diffeo- 
morphisms, described by a nonlinear sigma model, which includes higher derivative 
terms required to obtain the Fierz-Pauli mass term. 

Using our non-perturbative definition of massive gravity in de Sitter space, we 
discuss the appearance of an enhanced local symmetry and a null norm state in 
the linearized theory at the special value of the graviton mass M to the Hubble 
parameter H ratio [3H [351 [361 ESI HQ], which in general D dimensions occurs at 
M 2 = (D — 2)H 2 . In particular, by studying full non-perturbative equations of 
motion, we argue that there is no enhanced symmetry in the full nonlinear theory. 
We then argue that in the full nonlinear theory no null norm state is expected to 
arise at the aforesaid special value, which in turn suggests that no ghost might be 
present for M 2 < (D — 2)H 2 and the theory might be unitary for all values of M 
and H with no van Dam-Veltman-Zakharov (vDVZ) discontinuity [HI I42) 5 I . We 
argue that this is indeed the case by studying full nonlinear Hamiltonian for the 
relevant conformal and helicity-0 longitudinal modes. In particular, we argue that 
no negative norm state is present in the full nonlinear theory. 



4 For earlier and subsequent related works, see, e.g., |H El IH HI El 13 HI El HDl IH1 H21 H31 H31 Ull 

[161 El [HI H3 EQl EU [22] , and references therein. 

5 Absence of the vDVZ discontinuity in massive gravity in Minkowski space via gravitational 
Higgs mechanism was argued in [27] . 
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The rest of the paper is organized as follows. In Sections 2 and 3 we discuss the 
gravitational Higgs mechanism in the de Sitter background, which results in massive 
gravity in de Sitter space with the Fierz-Pauli mass term for the appropriately tuned 
cosmo logical constant. In Section 4 we discuss the enhanced local symmetry of the 
linearized theory and the absence thereof in the full nonlinear theory. In particular, 
we argue that the full nonlinear theory does not admit solutions obtained by trans- 
forming the de Sitter metric via such enhanced local symmetry transformations. In 
Section 5 we study the full nonlinear Hamiltonian for the relevant conformal and 
helicity-0 longitudinal modes and argue that no ghost is present for any values of 
the graviton mass and the Hubble parameter. We briefly summarize our conclusions 
in Section 6. 



2 De Sitter Solutions 

The goal of this section is to obtain massive gravity in de Sitter space via gravita- 
tional Higgs mechanism. Consider the induced metric for the scalar sector: 

Y M n = Z AB V M $ A V N $ B . (2) 

Here M = 0, . . . ,(D — 1) is a space-time index, and A — 0, . . . , (D — 1) is a global 
index. We will choose the scalar metric Z AB to be conformally flat de Sitter metric: 

Z AB = w 2 (0) tiab , (3) 

where 

*' s ss* • (4) 

v is a dimensionless coupling, and n A = (1, 0, . . . , 0). Also, let 

Y ee Y MN G MN . (5) 

The following action, albeit not the most generajj, will serve our purpose here: 

S Y = M°- 2 I d D xV^G [R - V(Y)] , (6) 

where a priori V(Y) is a generic function of Y. 



6 In the de Sitter as well as Minkowski cases one can consider a more general setup where 
the scalar action is constructed not just from Y, but from Ymn, Gmn and eM ...M D _ 1 , see, e.g., 
[26l Ell HOI El] . However, a simple action containing a scalar function V(Y) suffices to capture all 
qualitative features of gravitational Higgs mechanism. In particular, if this function is quadratic 
as in (PU|) the cosmological constant A must be negative in the context of Minkowski background 
(but not in the de Sitter case - see below), generically there is no restriction on A, which can be 
positive, negative or zero even in the context of the Minkowski background, once we allow cubic 
and/or higher order terms in V(Y), or consider non-polynomial V(Y). 



The equations of motion read: 

V M (V'(Y)Z AB V M <i> B ) = 9 -^- YV'(Y) , (7) 

Rmn — jGmnR — V kXT^mn — -^GMNViy) , (8) 

where prime denotes derivative w.r.t. Y. Multiplying ([7]) by Z^b^ s'P 3 an d con- 
tracting indices, we can rewrite the scalar equations of motion as follows: 



On, 



M 



V^GV'(Y)G MN Y NS ] - ]^GV'{Y)G MN d s Y MN = . (9) 



Since the theory possesses full diffeomorphism symmetry, (Q and (jHJ) are not all 
independent but linearly related due to Bianchi identities. Thus, multiplying (JS]) by 
\/—G, differentiating w.r.t. V N and contracting indices we arrive at ([9]). 
We are interested in finding solutions of the form: 

<\> A = m 5 A M x M , (10) 

G M n = exp(2A) t]mn , (11) 

where m is a mass-scale parameter, and A = A [umx m ) only depends on the time 
coordinate. The scalar equations of motion (|7|), or equivalently (jUj), then imply that 

A = - In (Hn M x M ) , (12) 

where H is the constant Hubble parameter. Furthermore, (jBJ) implies that 

Rmn — -GmnR = —-GmnA , (13) 

where 

A = V(YJ - 1 Y*V'(Y*) (14) 

is the cosmological constant for the de Sitter space whose metric is given by ffTTl) . 
and we have: 

y* = Dv 2 H 2 = D/2 2 , (15) 

H '' = (D - 1KB - 2) • (16 > 

Thus, the solution is completely determined by the "potential" V(Y) and the cou- 
pling v as Y"* is a solution of the following equation: 

Y * = (D - 1)(D - 2) [DV{Y * ) " 2Y * yl{Y ^ ■ W 

Note that, since A must be positive, the potential V(Y) cannot be completely arbi- 
trary. Also, Y"* and A are independent of m in ( ITUl) . which is due to the invariance 
of Y MN under simultaneous global rescalings 4> A — > \(fi A . 



3 Massive de Sitter Gravity 

In this section we study linearized fluctuations in the background given by (fTUj) and 
(iTTj) . Since diffeomorphisms are broken spontaneously, the equations of motion are 
invariant under the full diffeomorphism invariance. The scalar fluctuations <p A can 
therefore be gauged away using the diffeomorphisms: 



5<p A = V M A r = m 8 A M i 



M 



18) 



However, once we gauge away the scalars, diffeomorphisms can no longer be used to 
gauge away any of the graviton components Hmn defined as: 



G 



MN 



Gmn + h 



MN 



where 



Gmn = exp(2A) tjmn 



(19) 



(20) 



denotes the background de Sitter metric. Moreover, we will use the notation h = 

G MN h 



MN- 



After setting (p = 0, we have 

v 2 2 ~ 

Ymn — t ?r^ Vmn — ^ Gmn 



(nsx s ^ 2 
' mn<- 



V = YmnG mn = v 2 H 2 [D-h + 



Y* - fx 2 h 



(21) 
(22) 



where the ellipses stand for higher order terms in Hmn- 

Due to diffeomorphism invariance, the scalar equations of motion ([Jj) are related 
to (JHJ) via Bianchi identities. We will therefore focus on ([8]). Let us first rewrite it 
as follows: 



Rmn 

2 



-Gmn 



/' 



R-A = 

G M nV'(Y) - G M nV'(Y 



-Gmn [V{Y) - V%)] 



Linearizing the r.h.s. of this equation, we obtain: 



R 



1 



MN 



G 



MN 



R-A 



M 2 



GMNh — (h 



MN 



+ 



where 



M 2 = n 2 V'(Yj - 2/j,*V"(Y*) , 

CM 2 = 2/iV(i;) . 

This corresponds to adding a graviton mass term of the form 

- ^ [Ch MN h MN - h 2 ] 



(23) 



(24) 



(25) 
(26) 



(27) 



to the Einstein-Hilbert action with the cosmological constant A, and the Fierz-Pauli 
combination corresponds to taking £ = 1. This occurs for a special class of potentials 
with 

V'(Y.) = ~Y*V"(YJ . (28) 

Thus, as we see, we can obtain the Fierz-Pauli combination of the mass term for 
the graviton if we tune one combination of couplings. In fact, this tuning is nothing 
but the tuning of the cosmological constant - indeed, f|28|) relates the cosmological 
constant to higher derivative couplings. 
Thus, consider a simple example: 

V = A + Y + XY 2 . (29) 

The first term is the cosmological constant, the second term is the kinetic term for 
the scalars (which can always be normalized such that the corresponding coefficient 
is 1 by adjusting the coupling u), and the third term is a four- derivative term. We 
then have: 

Y - = ~2(3TT) A " • < 30 > 

which relates the mass parameter \x to the higher derivative coupling A: 



.2 



/' 



2(D + 2) 
and the graviton mass is given by: 



Y*/D = --— — - A" 1 , (31) 



M 2 = A" 1 . (32) 

(L> + 2) 2 v ; 



Note that we must have A < 0. Moreover, we have: 

D 2 + AD-8 , 



A = A t ^r- ^ ■ (33) 

4(£> + 2) 2 v ; 

Recall, however, that we have ( ITT)) . This implies that 

A- P-l)P-2) r i f34 , 

A " mdT2) a ' (34) 

and the cosmological constant A needs to be tuned against the higher derivative 
coupling A. 

Finally, note that 

A _ (D-l)(D-2)(D + 2) 

M2" 4^ • [6b) 

This ratio will become important in the next section. 



4 Enhanced Symmetry? 



In the previous section we saw that at the special value of the cosmological constant 
(given by (128]) ) we have massive gravity in de Sitter space with the Fierz-Pauli mass 
term in the linearized approximation. However, as we will argue in this section, 
there appears to be a qualitative difference between the linearized approximation 
and the full nonlinear theory, at least for some values of the coupling v. 

Before we do this, however, let us briefly comment on the counting of the prop- 
agating degrees of freedom. We started with massless gravity with D(D — 3)/2 
propagating degrees of freedom plus D scalars. However, just as in the case of mas- 
sive gravity in Minkowski space via gravitational Higgs mechanism discussed in [26J , 
at the spacial value of the cosmological constant ( T2"8"j) . due to the presence of higher 
derivative terms for the scalars, the kinetic term for scalar fluctuations reorganizes 
into that of a vector boson, and we have only D — 1 propagating scalar degrees of 
freedom. These D — 1 scalar degrees of freedom are eaten by the graviton in the 
process of spontaneous breaking of diffeomorphisms, the graviton acquires mass, and 
has (D + 1)(D - 2)/2 (= D(D - 3)/2 + (D - 1)) propagating degrees of freedom. 

Thus, in gravitational Higgs mechanism massive gravity arises as a result of 
spontaneous breaking of diffeomorphisms (as opposed to explicit breaking thereof 
by simply adding a mass term for the graviton fluctuations). In fact, gravitational 
Higgs mechanism provides a non-perturbative definition for massive gravity in the 
corresponding background. Once we gauge away the scalars via ([TBI , we obtain the 
following action for gravity in de Sitter background: 



S, 



G 



M°~ 2 



d D x\[^G 



R 



V(G MN G MN ) 



(36) 



where we have defined V(() = V(fi 2 (), and Gmn is the background de Sitter metric 
( II ip . The equations of motion are given by 

Rmn — -^GmnR = V{G Gkl)Gmn — t^GmnV{G ^ Gkl) , (37) 

and the Bianchi identities imply that 



d 



M 



V=G V\G KL G KL )G MN G 



NS 



V^G V\G KL G KL )G MN d s G M N = . (38 



Note that this condition is due to the presence of the "mass term" in (I3"S"|) . 

Let us now study linearized equations of motion. We expand Gmn = Gmn + 

Hmn'- 



□ h MN + VmV n H - V M V s h SN - V N V s h SM - G MN ph- V s V R h SR ] 



H 2 



2h<MN + (D — 3)GMNh 



-M A 



hiuN — GmnIt, 



(39) 



where Vm is the covariant derivative in the de Sitter background metric Gmn, and 

7 mV n . 



□ = G MJV Vj,fVftf. Furthermore, the condition (1551) reduces to 



V N h MN -V M h = . (40) 



Note that (HDD follows from 

While the linearized equations of motion (131?j) are not invariant under diffeo- 
morphisms, at the special value of the ratio A 2 /M 2 they are invariant under the 
following infinitesimal transformations: 

Sh MN = (Va/Vjv + H 2 G M n)x • (41) 

Indeed, (1391) are invariant under (14T1) when 

M 2 = (D- 2)H 2 . (42) 

The following identities are useful in deriving this result: 

(□ V M - V M □) X = {D - 1)H 2 V M % , (43) 

(D V M Viv - VmVjv □) X = 2H 2 (p VmV^v - G MN n) x • (44) 



In D = 4 the presence of the symmetry ( 14TT) at the point ( l4*2l) was discussed in [34J. 

The presence of this additional local symmetry in the linearized theory implies 
that at the point (1421) the graviton has (D + 1)(D — 2)/2 — 1 propagating degrees of 
freedom, one fewer than at generic points in the parameter space. Furthermore, the 
linearized theory is non-unitary for M 2 < (D — 2)H 2 (as the helicity-0 graviton mode 
has negative norm), while for M 2 > (D — 2)H 2 all (D + 1)(D — 2)/2 graviton modes 
are propagating and have positive norm [3U [351 HE] . Also, in the linearized theory 
at M 2 = (D — 2)H 2 the graviton Hmn can only couple to traceless conserved energy- 
momentum tensor Tmn as the coupling to the trace part of Tmn is inconsistent with 
the symmetry (I41[) . 

However, here we will argue that the additional local symmetry ( I4ip at the 
point (|42|) . which implies that the helicity-0 graviton mode has null norm at M 2 = 
(D — 2)H 2 and further acquires negative norm for M 2 < (D — 2)H 2 , is absent in 
the full nonlinear theory. To see this, let us start with the de Sitter solution and 
transform it via (14 ip with \ = at, where a is a constant, i.e., x is a linear function 
of time t = nsx s only. For such x, the transformation (1411) reads: 

This then implies that the so transformed metric is diagonal and of the form 

Goo = Goo, (46) 

G u = G lt {l + 2H 2 at) , (47) 



and is equivalent to the background de Sitter metric Gmn- We will now argue that 
the full non-perturbative equations of motion do not possess such a symmetry or 
such solutions. 

The following discussion can be straightforwardly generalized to general V(Y). 
However, for our purposes here it will suffice to consider quadratic V(Y) of the form 
( 129]) . Assuming ( jMJ) , from the previous section we then have: 



Au 2 = (D-2)(D + 2) 

A: 



D{D -A) , 

(X 



v(0 



2(D + 2) 



2{D + 2) 



[ jD(jD _4)_2( j D + 2)C + C 2 ] 



(48) 
(49) 

(50) 



Note that the cosmological constant A = in D = 4. 

For our purposes here it will suffice to first consider (138]) as opposed to the full 
equations of motion ( 137]) . We will look for solutions of the form 



G 



MN 



diag [G w J(t)G 



where /(£) is a function of time t only. For such solutions, we have 



V\G MN G 



MN i 



(D - 1)^ 



D + 2 

and (138]) reduces to the following equation for /(£): 
D+l 



f- 



D + l 
D- 1 



(51) 



(52) 



d t 



^('-§^)]+^(/-§^> + Ci>-iW-o. 



(53) 



Using the following equation 



d t V- 


-G = -V-G 


D + 


D-I3 t f 
2 / . 


•) 


we then have: 




i[(D + l)-{D- 


■3)f]d t f + ^^f(f-l) 


D + l' 
D-l. 


or equivalently: 




1 

2 


1 

_7" 


2 1 

f 1 ' f D+l 
J J D-l _ 


d t f + - t =o, 


with the solution given 


^y 













(f-iy 



f\f 



D+l I 
D-l I 

8 



(7^ 



(54) 



(55) 



(56) 



(57) 



where 7 is an integration constant. We therefore have (here we are looking for 
solutions where < / < (D + 1)/(D - 1)): 



./' 



D + l 

2(D-i 



iltf + eit 



D-l 



^m-^ 



I [1 + ( 7 t) 2 ] , (58) 



where e = ±1 and can be absorbed into the definition of 7: ej — > 7. Note that this 
solution is indeed bounded: < / < (D + l)/(D — 1). 

If we linearize (l5"6"j) via f — I + i/j, then we obtain the linearized solution (|47p . 
The solution (|58|) is therefore the non-perturbative counterpart of the perturbative 
solution ( |47|) with the identification 7 = — \J2(D — l)H 2 a. However, we will now 
show that the non-perturbative solution I J581) satisfies the full equations of motion 
(J37J) only for 7 = 0. 

To see this, it will suffice to consider the trace part of the equations of motion 
iB: 

1 \dv(0-2(v'(() , 



R 



D-2 



(59) 



where ( = G mn G M n = 1 + (D - l)f = D + (D - l)if>. Using J5QJ we then have: 



R 



H 



[(4 - D)C 2 + 2{D - 2){D + 2)C + D 2 (4 - D 
D(D-l)H 2 



D-l (D-1)U-D) o 



D 



8D 



(60) 



In D = 4 we have a simplification; however, we will continue to work in general D. 
From the definition of the Ricci scalar, we have: 



R = D{D-l)H 2 



i + a.h(/)-ifi?in(/) + l(a.iii(/)) 2 



(61) 



where 9* = tdt- We therefore have the following equation of motion for ip = f — 1: 



dJn(f)--d 2 Hf) + -(dJn(f)) 



D-l (D- 1)(4-D) 



D 



4 



D 



-0 



8£> 



v^ 2 . 



(62) 



Note that (162|) is exact. 

Note that the linearized version of (|62l) 



9*V - -^^ = — ^— ^ 



(63) 



admits solutions ^ = ^t, where £ is an integration constant, which match such linear 
solutions of the linearized version of (1581) . However, this does not hold beyond the 
linearized level. Indeed, it is not difficult to show that f[5"8"j) does not satisfy (jB^j) 



except for 7 = 0. A simple way to see this is to solve ( 162]) to the second order in t 
and compare this solution to (1581) expanded to the same order: 

Eq.fED: ^=^t+^-m 2 + Om 3 , (64) 

Eq.flBJ): ^ = V t _ (rjtf + 0( V t) 3 , (65) 

where we have defined 7 = \J(D — l)/2 7, and 77 is an integration constant. As we 
see, the two solutions do not match. Once we identify 77 = 7, the difference between 
the two solutions is given by: 

1* ~ 1>i = P 8 + p ( _^ 2) (^) 2 + °^f ■ ( 66 ) 



This implies that the full non-perturbative equations of motion ( 157|) do not admit 
solutions of the form floTj) (except for / = 1). 

An intuitive way of seeing why the linearized theory possesses the additional 
symmetry not present in the full nonlinear theory is as follows. In the linearized 
theory, (155]) reduces to fj4*0|) , which in turn implies that the part of the Einstein tensor 
containing covariant derivatives is traceless. The remaining part then is traceless 
at the special point (142]) . The tracelessness of the Einstein tensor is what leads to 
the appearance of the enhanced local symmetry in the linearized theory. Indeed, 
as we mentioned above, a conserved energy-momentum tensor for matter sources is 
compatible with this enhanced symmetry only if it is traceless. 

On the other hand, in the full nonlinear theory ( 155]) . unlike ( 140]) . depends on 
the structure of the potential V, and therefore it cannot possibly make the deriva- 
tive part of the Einstein tensor traceless as the latter intrinsically knows nothing 
about the structure of V. Therefore, non-perturbatively, one does not expect to 
have an enhanced local symmetry that would remove a propagating degree of free- 
dom. Indeed, for generic potentials V, the action (156]) does not possess any local 
symmetries. 

Let us quantify the previous two paragraphs. First, note that the scalars (fi A (or, 
their fluctuations around the background (fTUI) and ( ITT]) ) constitute the matter fields 
in the action ([6]). The energy-momentum tensor reads: 

T MN = -2M°~ 2 VGV'(Y)Y MN , (67) 

which is not traceless in the Higgs phase. Also, a nonlinear completion of (|4ip . 
before gauging away the scalars, is given by 

1 
SGmn = VmVatx + -zYmnX , (68) 

where Vm is the covariant derivative in the metric Gmn- in the Higgs phase, once 
we gauge away the scalars, we have 

8G MN = VmVjvX + H 2 G M nX , (69) 

10 



where Vm is the covariant derivative in the metric Gmn- However, ( 169|) is not a 
symmetry of the full nonlinear action ( 13 6 p . 

Here the following remark is in order. In ( 15 ip we assume that G 00 = G 00 . Here we 
can ask if there exist more general solutions with G 00 = g(t)G 00 (and G n = f(t)G n ) 

matching the perturbative solutions (J4U|) and ( HTp l. For this to be the case, we must 
have 

g = 1 + 0( V t) 2 . (70) 

However, if the symmetry ( 14 ip is indeed present, then we can always transform the 
metric such that g — 1 -+- 0(i]t) 3 . Indeed, for general x(t) that depends on time t 
only, (J4"T|) reads: 



5huN — ^M n N 



d tX + -dtX 



+ H 2 [td tX + x}G M N ■ (71) 



This implies that any 0(i]t) 2 term in g can be transformed away by including an 
appropriate 0(rjt) 2 term in \. More concretely, if g — 1 + a(r]t) 2 + 0(i]t) 3 , then the 
transformation (J4T|) with \ = arj 2 /3H 2 will result in g — 1 -+- 0(i]t) 3 . Therefore, if 
the enhanced local symmetry (|4T|) were present, our assumption of G 00 = G 00 would 
hold. 



5 Is there a Ghost? 

In the previous section we argued that in the full nonlinear theory there is no 
enhanced local symmetry at M 2 = (D — 2)H 2 . Absent such local symmetry, one does 
not expect appearance of a null norm state at M 2 = (D — 2)H 2 , which would make 
it difficult to imagine that a negative norm state would appear for M 2 < (D — 2)H 2 . 
This suggests that no ghost might be present for lower graviton mass values and 
the theory might be unitary for all values of the graviton mass M and the Hubble 
parameter H with no vDVZ discontinuity^] The purpose of this section is to argue 
that this is indeed the case. We will do this by studying the full nonlinear action 
for the relevant modes, which we identify next. In particular, we will argue that no 
negative norm state is present for any value of the ratio M/H. 

To identify the relevant modes in the full nonlinear theory, let us note that in the 
linearized theory the potentially "troublesome" mode is the longitudinal helicity-0 
mode p. However, we must also include the conformal mode u as there is kinetic 



7 Based on symmetry considerations, namely, the SO(D— 1) invariance in the spatial directions, 
the off-diagonal terms in G MN are not relevant in this discussion. 

8 The related issue of causality will be relegated to future work. The possibility of violating 
causality by means of superluminal propagation of signals, which affects other models of modi- 
fied gravity, cannot be checked easily in the model at hand. Namely, the phase velocity at low 
wavelengths of a linearized perturbation, the speed of propagation of a signal, has no meaning in 
a model where, as argued above, the linearized analysis is not applicable, and the study of the 
corresponding phenomenon at nonlinear level is beyond the scope of this article. 

11 



mixing between p and u. In fact, p and u are not independent but are related via 
Bianchi identities. Therefore, in the linearized language one must look at the modes 
of the form 

h M N = G M n u) + V M V N p . (72) 

Furthermore, based on symmetry considerations, namely, the SO(D — 1) invariance 
in the spatial directions, we can focus on field configurations independent of spatial 
coordinates. Indeed, for our purposes here we can compactify the spatial coordinates 
on a torus T D_1 and disregard the Kaluza-Klein modes. This way we reduce the D- 
dimensional theory to a classical mechanical system, which suffices for our purposes 
here. Indeed, with proper care (see below), if there is a negative norm state in the 
uncompactified theory, it will be visible in its compactified version, and vice-versa. 
Let us therefore consider field configurations of the form: 



G 



MN 



dmg(g(t)G m J(t)G 



where g(t) and f(t) are functions of time t only. The action 
follows: 



S, 



G 



— K 



^g-h-^{W + V(g + n) 



(73) 
then reduces as 

(74) 



where 



Q = (D-l)f , 



(75) 

(76) 

(77) 



and Wd-i is the volume in the spatial dimensions (i.e., the volume of T D_1 ). Note 
that g is a Lagrange multiplier. The goal is to integrate out g and obtain the 
corresponding action for /. It is then this action that we should test for the presence 
of a negative norm state. 

The equation of motion for g reads: 



V(g + ft) - 2gV'(g + ft) = AgU 2 . 



(78) 



The following discussion can be straightforwardly generalized to general V. However, 
for our purposes here it will suffice to consider quadratic V corresponding to (129]) . 
We then have: 



3Xp 2 g 2 



1 + 2\p 2 tt + —U 
p 2 



g -^ + n[i + \p 2 n]\ =0. 



(79) 
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We can therefore express g in terms of / and <9* ln(/): 

A 



6A/x g 



i + 2\fi z n + — u z 
/i 2 



+ 



u 



I + 2\fl 2 tt + —U 2 

n 2 



-i 2 



12A/i 2 { 4 + [1 + A/i 2 fi] 
/i 2 



(80) 



where the branch is fixed by the requirement that g = 1 when / = 1. Substituting 
the so expressed g into (171)1 . we obtain an action which is a nonlinear functional of 
/ and <9*ln(/). 

For our purposes here it is more convenient to work with the logarithmic time 
coordinate r and the canonical variable q, where 



r = \n(Ht) , 
q = \n(f) + 2t , 
Q = {D- l)e q ~ 2T 

U = -d T q , 



and the action reads: 
So — I drL 



■kH 



''-' I dr g-h-^ 



AgU 2 + V(g + Q) 



(81) 
(82) 
(83) 

(84) 



IB) 



where L is the Lagrangian. This action corresponds to a classical mechanical system 
with a lagrange multiplier g and a time-dependent potential. Upon integrating out 
the Lagrange multiplier, the time dependence also propagates into the "kinetic" (or, 
more precisely, momentum-dependent) terms. 
Next, the conjugate momentum is given by 



dL 



P 



d{d T q) 



-kH 



D-l. 



X 



< < ^g^gAU 2 - )^g-igV(g + Q) + g-*gV'(g + ft)+ g^AU 



where 
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dg 



d(d T q) ' 
Using ([78]), (PJ simplifies to 

p = -nH^e-^g^AU 
and the Hamiltonian is given by 



(86) 
(87) 



H=pd T q- L = -KH u ~ L g-*e 



AgU 2 -V(g + Q) 



19) 
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We can now see if this Hamiltonian is bounded from below. 
First, using ( 175|) . we have: 



D-l„l„-i 



U = 2KH n - 1 g2 e -^ q V'(g + O) = 2^KH u - l g*e 



[l + 2\fi 2 (g + n)] . (90) 



Using 



where 



we can rewrite this Hamiltonian as follows: 

H = -^KH^gh-^ [X - Z] 



(91) 



X 



A 



P 



1 + 2A/i 2 fi + — U 2 
[i 1 



U 2 - 4A/z 2 ft - 2 . 



-I 2 



- i2A// 2 <j4 + ^[H-A/i 2 fi]^ ■ (92) 



(93) 



The presence of a ghost would imply that the Hamiltonian is unbounded from below 
for large values of U 2 (recall that U 2 contains the "kinetic" term). However, it is 
not difficult to show that this Hamiltonian suffers from no such pathology. Indeed, 
we can rewrite it as follows: 



n 



~pL 2 KH D -Y*e 



1 D-l 



X 2 



x + z 



2^KH u ~'g^e 



2j_i i _£=!« 4AA - 1 + 2(A//i 2 ) U 2 [1 + 2A/i 2 fi] 



X + Z 
which in the large U 2 limit reads: 

H = 2 f i 2 KH D ~ 1 gh- I ^ lq [1 + 2X/j 2 n] + 0(1/U 2 ) . 

Furthermore, from ( 180]) we have 



6Xfi 2 g 



X 2 -Q 2 
X + Q 
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AA + Xfi 2 Q [1 + A/i 2 ft] 
X + Q 



0{l/U 2 



where 



A 

7 2 



Q = ^U 2 + 2Xfi 2 Q + 1 . 



(94) 



(95) 



(96) 



(97) 



So, in the large U 2 limit the Hamiltonian actually vanishes. 

Note that the above argument implicitly assumes that Q is bounded from above. 
This is indeed the case as g must be at least non-negative, which implies that 



AA + X/j, 2 n [I + Xn 2 tt] < 



(9* 
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and Q is bounded as follows (note that we must have £1 > 0): 



1 - Vl-4AA\ n 1 + y/l - 4AA . . 



and we must further have 



A > L , (100) 



which together with ( 155]) implies that we must have 

1 £ (flf^ 1 • (101 > 

which is always satisfied as A < 0. 

Thus, as we see, there appears to be no ghost in the full nonlinear theory. To 
understand why a ghost is present in the linearized theory, let us linearize our 
Hamiltonian. To do this, we will assume that / = 1 + ip, where \ip\ <C 1, i.e., we 
consider small fluctuations around the de Sitter background. Furthermore, we also 
assume that \d T ip\ C 1, so we can linearize the Hamiltonian to the second order in 
d T i/) as well as ip. In fact, here we are interested in the terms containing (d T ip) 2 . A 
straightforward computation gives the following linearized Hamiltonian (the ellipses 
stand for the terms not containing (d T ip) 2 ): 

H* = -^ K H D - l e- {D -^ T (d T ijf + ... , (102) 

i.e., this linearized Hamiltonian contains a ghost for all values of M and H, which 
is absent in the full nonlinear theory. In this regard, we discuss a simple illustrative 
example in Appendix lAl 

Here one might find it puzzling that, all the potential pitfalls of linearization 
notwithstanding, in the linearized theory the ghost appears for all values of M and 
H, while according to [5U [351 EH] a negative norm state is expected to appear 
only for M 2 < (D — 2)H 2 . The difference here is due to the parametrization of 
the conformal and helicity-0 longitudinal modes. We have been working with (175]) . 
while the aforesaid result of [M] |5J2 155] applies to (172]) . The difference between the 
two is that ([75]) has no derivatives. In this regard, one might wonder if the ghost 
is "masked" by (IT5|r l. Furthermore, one might wonder if the dimensionally reduced 
action ([85]) provides an adequate description. 



9 Considcr a ghost in D dimensions: L = —c d M 4>dM<j>, c < 0. Compactify on T D ~ X : L = 
ci(d t (j)) 2 (c\ = cWd-x, Wd—i being the volume of T ^ 1 ). The Hamiltonian is not bounded from 
below: % — p 2 /4ci, where the conjugate momentum p = 2c\dt4>. However, let q = dt<j>. If q 
is treated naively as the canonical variable, then L = c\q 2 , "H = —c\q 2 , which is bounded from 
below. That is, the ghost appears to have been "masked" by redefining the variables. The flaw 
in the argument is that the transformation of variables is not a canonical one, therefore, the new 
Hamiltonian describes a different dynamical system. 
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In this regard, we have explicitly checked that if we expand the action ( 1851) to the 
quadratic order in the parametrization corresponding to (1721) . we obtain that there 
is a null norm state at M 2 = (D - 2)H 2 and a ghost at M 2 < (D - 2)H 2 , so the 
dimensionally reduced action f )85|) correctly reproduces the linearized results [37] in 
the parametrization corresponding to (T72l) . In fact, what transpires is the following. 
At the quadratic order, the second derivatives introduced by the parametrization 
corresponding to (1721) can be integrated by parts to arrive at an action containing 
only first derivatives of u and p. This action then possesses the aforesaid prop- 
erties w.r.t. the appearance of a null norm state and a ghost. However, we have 
explicitly checked that already at the cubic level the second derivatives introduced 
by the parametrization corresponding to (172]) cannot be integrated by parts, so 
the resulting action invariably includes terms with second derivatives of p. This is 
clearly problematic already at the cubic level and suggests that the parametrization 
corresponding to (172]) cannot be used beyond the linearized approximation. 

Indeed, according to (1721) the ghost would appear for M 2 < (D — 2)H 2 , while at 
M 2 = (D — 2)H 2 we would have a null norm state, and for M 2 > (D — 2)H 2 the 
theory is unitary. The appearance of a null norm state would signal the presence 
of an enhanced local symmetry at M 2 = (D — 2)H 2 . However, in Section 4 we saw 
that there is no such enhanced symmetry in the full nonlinear theory. Therefore, in 
the full nonlinear theory we either have a ghost for all values of M and H, or the 
theory is unitary 10 ! for all values of M and H. This suggests that the "special" point 
M 2 = (D — 2)H 2 arises in the linearized theory both due to linearization and the 
parametrization (1721) . Indeed, a nonlinear completion of ( 1721) is given by: 

Q MN = qMNj + V M V N U (1Q3) 

where the covariant derivative is defined w.r.t. the metric Gmn (this choice does not 
affect our discussion here). Note, however, that such a parametrization of the metric 
is rather problematic in the context of the full nonlinear theory as it introduces higher 
derivative terms in u, which should therefore not be used as the canonical variable 
in the full nonlinear theory. This suggests that our parametrization (1731) is indeed 
adequate 11 !. As we saw, in this parametrization there is a ghost in the linearized 
theory for alO values of M and H. However, the full nonlinear theory appears to 
be unitary for all values of M and H. 

Before concluding, let us comment on another related issue. In the linearized 
theory, a "fifth constraint" (in the D = 4 language), which is complementary to 
the Bianchi identities, removes one degree of freedom (see, for example, [3S] for a 



10 In this regard, note that, on general grounds, for M 2 3> H 2 one expects no ghost to be present. 

n In fact, without giving details, let us simply mention that, one arrives at the same conclu- 
sion by adding a Lagrange multiplier r\ leading to a constraint g = / — H 2 d 2 u, which is a non- 
linear completion of (|72p . and which follows from the following additional term in the action: 
-kH - 1 J r { V [g - f] - H 2 d TV d T u}. 

12 In this regard, there is no ghost "masking" as a ghost does arise in this parametrization upon 
linearization. 
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derivation). At the special point M 2 = {D — 2)H 2 , this constraint trivializes and 
leads to an enhanced symmetry which removes yet another degree of freedom - as 
mentioned above, the norm of one of the modes becomes null. At the nonlinear level, 
however, these features disappear since extra momentum and position dependent 
terms also contribute to the would-be constraint. In this regard, one might worry 
whether a propagating (ghostlike) "sixth mode" (in the D = 4 language) might also 
be present. However, our non-perturbative Hamiltonian analysis appears to indicate 
that no ghostlike states or instabilities are present in the full nonlinear theory. 

6 The Upshot 

As we argued in Section 4, in the full nonlinear theory there is no enhanced local 
symmetry at the special point ( 14*2|) . In particular, we argued that the full nonlinear 
theory does not admit solutions obtained by transforming the de Sitter metric via 
such enhanced local symmetry transformations. 

Absent such enhanced local symmetry at the special point (1421) in the full nonlin- 
ear theory, we do not expect to have a null norm state and a reduction in the number 
of propagating degrees of freedom. Furthermore, there is also no reason to believe 
that for M 2 < (D — 2)H 2 there is a negative norm state in the full nonlinear theory. 
If so, we can expect that, in the context of gravitational Higgs mechanism, where 
diffeomorphisms are broken spontaneously, there should be no van Dam-Veltman- 
Zakharov (HI H2] discontinuity, absence of which in the context of massive gravity 
via gravitational Higgs mechanism in Minkowski space was argued in [27]. In fact, 
from our analysis of the full nonlinear Hamiltonian for the relevant conformal and 
helicity-0 longitudinal modes it indeed appears that there is no ghost in the full 
nonlinear theory. In this regard, it appears that the vacuum corresponding to the 
linearized theory is just another vacuum that is unstable, and that the true vac- 
uum corresponding to the full nonlinear theory is a different and, apparently, stable 
vacuum. Here one analogy that comes to mind is ghost condensation: An approxi- 
mated version of a theory appears to have a ghost, but in the true vacuum the ghost 
condenses and the theory lacks such pathologies. 

The upshot is that we have presented evidence that gravitational Higgs mecha- 
nism may provide a non-perturbative definition of massive gravity (both in Minkowski 
and de Sitter backgrounds, as well as in general curved backgrounds [43J ) , and it ap- 
pears that massive gravity obtained via gravitational Higgs mechanism, on general 
grounds, may be consistent since diffeomorphisms in gravitational Higgs mechanism 
are broken spontaneously (as opposed to explicit breaking by simply adding the 
Fierz-Pauli term for the graviton). This may open a new arena for studying infrared 
modified (that is, massive) gravity in the context of cosmology with non-vanishing 
cosmological constant. 
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A Some Pitfalls of Linearization 

In this appendix we illustrate some pitfalls of linearization. Namely, apparently a 
nonlinear theory can have a bounded-from-below Hamiltonian, while a linearized 
version thereof can have a ghost. Consider the following simple toy Lagrangian for 
a single scalar field (f>: 

L = cVY 2 + 2aY + b 2 , (104) 

where 

Y = d M <pd M <P , (105) 

b > 0, and a can be either positive, negative or zero. To make sure that the square 
root in f)104|) is well defined, let us assume that 

\a\ < b . (106) 

Next, suppose we naively linearize (to the quadratic order in or, equivalently, the 
first order in Y). The so linearized Lagrangian reads: 

U = cb+(ca/b)Y + 0(Y 2 ) , (107) 

which has a ghost for 

ca > . (108) 

However, the full nonlinear theory does not possess a ghost. To see this, let us first 
compactify the spatial directions on a (D — l)-dimensional torus T D ~ X and reduce 
the theory to a classical mechanical system. The Lagrangian now reads: 



L = Cl VZ 2 -2aZ + b 2 , (109) 

where c\ = cWd-i, Wd-\ is the volume of T D_1 , and 

Z = {d t <pf (110) 

The conjugate momentum is given by 

,- y *-">(*»> , (m) 

VZ 2 - 2aZ + b 2 
and the Hamiltonian reads: 

H= »<*•-*> . (112) 

^Z 2 - 2aZ + b 2 

Because of (11Q6J) . the square root is well defined. Furthermore, for c > (i.e., ci > 0) 
the Hamiltonian is bounded from below regardless of the sign of a. Therefore, there 
is no ghost in the full nonlinear theory, and its appearance in the linearized version 
thereof is due to linearization. 
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